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ABSTRACT 

Vertical plane response of submersible vehicles in the proximity of a free surface in 
deep water is evaluated using a potential flow, strip theory solver. Two criteria, that are 
periscope submergence, and sail broaching are used to quantify the response. These criteria 
combined with the vehicle’s response amplitude operators in regular sinusoidal waves along 
with a statistical description of the seaway lead to an assessment of an overall operability 
index for the vehicle. This thesis presents a systematic parametric study of the effects of body 
geometry on near surface response. Two cases, namely limited diameter and limited length are 
considered. The total volume of the vehicle is kept constant, and certain shape factors are 
changed, while either the overall diameter or the overall length remains the same. The 
operability index is calculated for each case within a given range for sea states and sea 
directions and for various shape factors, vehicle speeds and operating depths. The results 
indicate that certain changes of shape factors can improve vehicle operations in various depth 
and speed combinations. 
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I. INTRODUCTION 



Submarine periscope depth operations are conducted in order to accomplish 
specific tasks and perform them in minimum time and with the least amount of noise. 
Accurate submarine maneuvering predictions are essential both for design to assess 
alternatives based on their maneuvering performance, and for operation, to provide 
optimal and safe submerged operating envelopes. Submarines are subject to several 
exciting forces and moments at periscope depth beneath a seaway. These forces and 
moments induce an oscillatory motion at the wave frequencies and second-order drifting 
motions at very low frequencies well outside the wave spectrum of the seaway, which are 
referred to as free surface suction effects. It is difficult in practice to control these low 
frequency motions and may result in a rise to unsatisfactory depth keeping. Standard ways 
of computing free surface effects rely on combinations of potential flow and semi- 
empirical coefficient based models. Rankine type sources are distributed along the hull of 
the ship, which satisfy the free surface boundary condition. Source strength can be 
computed by satisfying the exact body boundary condition, that no fluid can pass through 
the hull surface. Discretization of the hull form into a finite set of Hess-Smith type 
quadrilateral panels allows formulation of algebraic system of equations to be solved for 
the unknown singularity strengths [Ref. 1]. Combination of forces and moments generated 
in this way with deep water force predictions can then be utilized to simulate the motion of 
the boat under waves. These forces and moments are slowly varying in time and as a result 
they can be controlled by either the operators or automatic control systems. In this work a 
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potential flow, strip theory solver program is utilized to determine vehicle motions in the 
proximity of a free surface in deep water, based on the work by Beck and Troesch [Ref. 
2 ]- 

Tactical assessment is possible by adopting a number of criteria, each pertaining to 
different operational hazards. These criteria can be divided into two main categories, such 
as subtle and catastrophic failures. In the study, we have considered two criteria which are 
periscope submergence and sail broaching. Periscope submergence is a subtle failure and 
sail broaching is a catastrophic failure. 

Subtle failures are the events which will occur in all types of periscope depth 
operations, such as, propeller emergence, mast emergence, periscope submergence, mast 
submergence. The frequency of these events imposes operability limits for a certain sea- 
state, although single occurrence of them does not constitute failure of operations. 
Periscope submergence impairs visual information. The dominant criterion is the number 
of occurrences per unit time. 

Catastrophic failures are the ones which will probably result in either cancel of 
operations, failure to complete mission or submarine detection, such as sail broaching, loss 
of depth control deep. Since a single event occurrence is identified as a catastrophic 
failure, the primary analytical criterion is the expected duration between two consecutive 
events, instead of being the number of events per time. Broaching is defined as the loss of 
depth control shallow to a sufficient extent and duration. It is assumed that the submarine 
detection will occur with probability of one each time a broaching occurs. 
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In Chapter II of the thesis, evaluation of the response is described, beginning with 
the definition of submarine geometry. A typical hull consists of three sections, the entrance 
(bow) the shape of a parabolic body of revolution, the parallel middle body of a cylindrical 
shape, and the run (stem) the shape of an ellipsoid of revolution (Jackson, 1992) [Ref. 3]. 
The parameters and all the coefficients are explained and formulated. Chapter III shows 
the resulting operability indexes of a submarine for different sea-states and wave heading 
angles and deep water operations for periscope submergence, sail broaching and 
combined criteria. These results are obtained from the limited diameter and limited length 
cases, different shape factors and speed/operating depth combinations. Conclusions from 
the study and recommendations for further studies are presented in Chapter IV. 
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II. EVALUATION OF RESPONSE 



A. GEOMETRY OF A SUBMARINE 

A typical submarine hull is considered as a generic body of revolution, which is 
rotated about a line parallel to the center-line, as described by Jackson (1992) [Ref. 3]. It 
has a length/diameter ( UD ) ratio of six and a maximum diameter at 0.4 L. The body is 
composed of three main sections, as shown in Figure 1 [Ref. 3]. The forward section is 
called the entrance, which is a portion of an ellipsoid of revolution. The middle section is 
the parallel middle body (PMB) with a cylindrical shape. The third section is the after end 
called the run, which is composed of a paraboloid of revolution. The entrance has a length, 
Lf , of 2.4 diameters. The run has a length, L a , of 3.6 diameters. The algebraic sum of the 
lengths, Lf, L a , and the length of the PMB , Lpmb , is the overall length of the hull. 




Figure 1. Submarine Geometry [Ref. 3] 
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The body coordinates, which define the forward and aft shapes are given by. 




( 1 ) 



D 




( 2 ) 



The values of Xj and x a are the offsets from the maximum diameter, and y f and y a are the 
radii at the respective offset points. 

The exponents ( n a and n/) in Equations (1) and (2) are the shape factor 
coefficients, which control the shape of the fore and aft bodies, respectively. Higher values 
of these coefficients correspond to fuller hull shapes, and lower values to finer shapes. The 
effects of changing the shape factor coefficients on the hull shape are shown in Figures 2 
and 3, where the hull shapes for three values of the shape factors 2, 3 and 4 are shown. In 
this study, it is assumed that the total volume of the ship remains the same (so that ship 
displacement does not vary), while either overall length or diameter remains the same. We 
refer to the first case as the limited length case or inactive diameter constraint, and to the 
second case as the limited diameter case or inactive length constraint. Both of these two 
cases will be analyzed in our parametric studies. 

If one were to use equations for true ellipsoids and parabolas, the entrance and the 
run would be too fine for a modem submarine. The displacement can be increased by 
using larger shape factors or higher Lpmb ■ The prismatic coefficients , C pa and C p j ,are used 
to calculate volumes. For a cylinder its prismatic coefficient is 1. For a submarine-like 
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Figure 2. The effect of changing the shape factors for the limited diameter case. 




Figure 3. The effect of changing the shape factors for the limited length case. 
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body the prismatic coefficient can be evaluated in terms of its geometry. Using the above 
concept, Lpmb is the difference between the overall length and six times the diameter, that 
is L-6D. A method of calculating the volume of the entire hull can be developed by 
calculating the volume of each section separately by using the expressions above. Let Vf, 
V a , and Vpmb denote volumes of the entrance, the run, and the PMB, respectively. The 
resulting equations are: 

v , =^(C„ 2 - 4£> ) (3) 

V,=^-(C^.6D) (4) 



^ PMB ~ 



jzD 



(L-6D) 



The above can be combined into the following. 



V = 



kD 



3.6C P a H 6 + 2.4G/ 

D 



(5) 



( 6 ) 



It can be easily shown that the prismatic coefficients can be calculated as. 



i 



C pf = J (1 - x"' ) 2/n/ dx 
0 


(7) 


1 

C pa =\(\-x n ‘) 2 dx 


(8) 



o 



These integrals are numerically evaluated using the built in “quad” function in Matlab, 
although analytic evaluation is possible. 
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(9) 



C, 



pa l + 3n a +2n 2 a 




( 10 ) 



where the Gamma function is defined in Ref. 4 (Abramowitz and Stegun, 1970). 

Initially, the total volume is calculated for the values of L =360 ft., D = 30 ft., 
n a = 3.0 and n/= 3.0, by using Equations (6) through (8), and yielded a value of 
217337.73 ft 3 . This value is kept constant throughout the calculations. 

For the limited length case, L is kept constant at 360 ft. The shape factors n a and 
rij are varied between 2.0, 3.0 and 4.0. The prismatic coefficients C pa and C p ; are 
calculated for each of them by evaluating the integrals in Equations (7) and (8) 
numerically. Then the corresponding diameters are found by solving for the maximum hull 
diameter D in Equation (6). This is achieved by solving the following cubic equation, 



Two iterations ensure that the solution converges to the value which meets the 
requirement to have L a =3.6D and Lf=2AD. 

For the limited diameter case, D is kept constant at 30 ft. C pa and C p j are 
calculated the same way as the limited length case. Then the corresponding lengths are 
found by solving for L in Equation (6).The resulting equation becomes. 



(3.6C P a + 2.4 Cp/ - 6)D 3 + LD 2 



4V 



= 0 



(ID 



K 




( 12 ) 
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Two Matlab programs (Appendix) are used to perform the calculations for each 
case. In the limited length case, the program “limlen” starts by inputting n a and «/, and 
computes the corresponding diameter. Similarly, for the limited diameter case, the 
program “limdia” computes the length for the given n a and n/ values. These values are 
used as an input in the strip theory seakeeping prediction program. 



B. MOTIONS IN A SEAWAY 

Wave patterns in an open sea are ever changing with time and space, in a manner 
that appears to defy analysis be it linear or second order Stokes [Ref 5]. Ambient waves 
on the surface of the sea are dispersive as well as random. Random refers to the character 
of the wave height distribution. In a continuous distribution, the sinusoidal waves have 
continuously distributed amplitude and phase so that in summation the variation of wave 
height with time is not systematic in any respect, but random. The generating mechanism 
is, predominantly, the effect upon the water surface of wind in the atmosphere. The 
practically useful data extractable from a random wave record h(t) is its spectral density, 
5(co). The random h(t) record is processed in such a way to produce a curve of S( (0) 
versus wave frequency, co. The spectral density is obtained from a wave height record 
taken over a time period for which the sea conditions are assumed to be unchanging, in an 
average sense (stationary). This corresponds to a certain sea state. The function S((o,Q ) is 
called the spectral energy density or simply the energy spectrum. More specifically, this is 
a directional energy spectrum; it can be integrated over all wave directions to give the 
frequency spectrum. 
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( 13 ) 



S(co) = S((0,d)de . 



Usually in the fields of ocean engineering and naval architecture it is customary to 
assume that the waves are long crested which means the fluid motion is two dimensional 
and the wave crests are parallel. With such a simplification it is possible to use existing 
information for the frequency spectrum (13), which is based on a combination of theory 
and full scale observations. The sea spectrum (spectral density) gives us information on 
mean wave height within finite frequency bands. Since most of the wave energy is within a 
relatively small range of wave lengths where it may resonate the ship, we can model the 
seaway as a narrow band random process. 

For most purposes we are interested primarily in the larger waves. The most 
common parameter that takes this into account is the significant wave height, H m , defined 
as the average of the highest one third of all waves. This is computed by 



In this equation, mo is the area under the spectrum 5((o) integrated over the entire range of 
frequencies CO. An average frequency of the spectrum can be defined as the expected 
number of zero upcrossings per unit time, that is, the number of times the wave amplitude 
passes through zero with positive slope. The final result here is 



i/,/3 =4.0 (m 0 )' /2 . 



( 14 ) 




V /2 



( 15 ) 



The average period between zero upcrossings is 




( 16 ) 
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More meaningful frequency parameters can be obtained from the set of moments, 



which depend on spectrum shape 



m n = j CO n S(coi)dco , n= 0,1,2,... 



(17) 



o 



In particular, the area, m 0 , is the variance or the total energy of the spectrum. Also m is 
variance of velocity and m A is variance of acceleration. 

A good model for fully developed seas is the classical Pierson-Moskowitz 
spectrum. This spectral form depends upon a single parameter which is the significant 
wave height. It is intended to represent point spectrum of a fully-developed sea. Fetch 
and duration are assumed to sufficiently large so that the sea has reached steady state, in a 
statistical sense. This spectral family should be recognized as an asymptotic form, reached 
after an extended period of steady wind, with no contamination from an underlying swell. 
Using the spectral family, along with the similarity theory of S. A. Kitaigorodskii, Pierson 
and Moskowitz (1964) [Ref. 5] arrived at the following analytical formulation for ideal sea 
spectra, 




\ 2 

8 



(18) 




where 



S, + (co) = one-sided incident wave spectrum 



g = acceleration of gravity 



H m = significant wave height 



0) = wave frequency 
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In Figure 4 we can observe typical Pierson-Moskowitz wave spectra for 5 m. significant 
wave height. 




Any conclusions drawn on the seakeeping behavior of a ship based on the critical 
examination of motion response in regular waves can, at best, assume only academic 
significance. The establishment of the seakeeping behavior of a ship has to be done in a 
realistic seaway. With the spectral description of sea waves given before, we can return to 
the subject of body motions and generalize the results of regular harmonic waves. If the 
sea waves are described by the random distribution , and if the response of the body to 
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each component wave is defined by a response amplitude operator Z(co,Q), the body 
response will be 

= . ( 19 ) 

The principal assumption here is that linear superposition applies, as it must in any event 
for the underlying development of the RAO and the spectrum. 

Like the waves themselves, the response (19) is a random variable. The statistic of 
the body response are identical to the wave statistics, except that the wave energy 
spectrum S is multiplied by the square of the RAO (this is a property of linear systems). 
Thus, if the subscript R represents any body response, we have 

S R {co) = \Z R {cof S(co) , ( 20 ) 

where Z„( on) is the RAO of the response R, and S(co) the spectrum of the seaway. Equation 
(20) can then be utilized to obtain the spectrum of the response R. Figure 5 displays the 
spectrum of response of the relative vertical motion at the top of our model submarine’s 
sail while submarine’s forward speed is 5 Knots and it is at 3 submarine diameter depth. 
Also seaway is modeled by Pierson-Moskowitz spectrum with 5 m. significant wave 
height and head seas. 

To a large extent, equation (20) provides the justification for studying regular 
wave responses. The transfer function Z*( to) is valid not only in regular waves, where it 
has been derived, but also in a superposition of regular waves, and ultimately in a 
spectrum of random waves. Generally speaking, a vessel with favorable response 
characteristics in regular waves will be good in irregular waves, and vice versa. 
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wave frequency (rad) 

Figure 5. Spectrum of response for relative vertical motion. 

The average period between zero upcrossings was determined by Equation (16), 
and the number between zero upcrossings per unit time is 



z 2ft 



( 21 ) 



where m * , m* are the moments of the particular response R, whose spectral density is 
given by Equation (20). Equation (21) can be generalized for the case of the average 
number of upcrossings above a specified level a as in 



Ka = T~lP^ eX P 

2 K\m 0 



a 



2 



2 m o j 



( 22 ) 
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Equation (22) can be utilized to determine such events deck wetness and bow slamming 
for a surface ship or periscope submergence and sail broaching for a near surface 
submarine. If/ represents height of the periscope over calm sea surface level, the number 
of periscope submergence events per hour is 



N = 3600— —exp 
p 2 k V ra n F 









2m, 



( 23 ) 



oj 



where m 0 , m 2 are the moments of the vertical relative motion spectrum at periscope. The 
same equation can be used to estimate the frequency of sail broaching, with / substituted 
by the distance between top of the sail and encountered wave surface. Of course, m„ , m 2 
are now the moments of the relative motion spectrum at top of the sail. 
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III. RESULTS 



A. OPERABILITY INDEX 

In this study two criteria are considered for the near surface operations of the 
submarine in deep water, for both the limited diameter and the limited length cases. First is 
the number of periscope submergence events per hour (N pI ), which can be calculated using 
Equation (23) and selected as 300, which corresponds to five periscope submergence 
events per minute. This is an arbitrary number and different numbers could be picked for 
different operational considerations. The same number is used for all cases to represent all 
possible choices. The other criterion is the number of sail broaching events per hour (N p2 ), 
which can be also calculated from Equation (23) and selected as one. Since a single event 
occurrence may result in submarine detection, this places far greater emphasis on sail 
broaching than periscope submergence. 

Having the tools to compute the two performance indices defined above in a given 
seaway. Suppose that the submarine conducts periscope depth operation in a seaway 
characterized by a significant wave height (Pierson-Moskowitz Wave Spectrum), so that 
the sea spectrum is defined. For all round the clock boat headings relative to the 
predominant wave direction for which the operations are to be conducted, a polar plot 
diagram similar to the one in Figure 6 is prepared. Significant wave heights are 
represented along the radial direction of the polar plot. The shaded area in the plot shows 
wave height and wave direction combinations where the selected tactical assessment 
criterion is exceeded. Three plots are combined on a single plot for different shape factors, 
which are selected as two, three and four at the same speed/depth combination. Letting 
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the polar area of the disk in Figure 6 be A„ and the subset of A„ within which the boat can 
conduct the operation be A, a performance index characterizing the ability of the boat 
perform this operation in the specified submarine velocity and depth can be defined as 
100(A/A„). Generally a submarine’s forward speed ranges from three to twelve knots in 
periscope depth operations. We used three, five, eight and eleven knots submarine forward 
speeds, U, in our calculations. Depths, h, beneath the surface were selected from 1.5 to 
2.5 boat diameters measured from the keel up. The calculations are repeated for both the 
limited diameter and the limited length cases and for the shape factors n a and n/, which are 
selected as two, three and four. In computing the above index we could easily take into 
account the probability of occurrence of a particular sea state and wave heading angle in 
the area of interest by introducing appropriate weight factors. In this study we assumed 
that all possible sea-states and wave heading angles are equally probable. In the following 
sections we discuss the results for both criteria/operability indexes and also for the 
combined criterion/operability index, where both criteria are taken into consideration at 
the same time. In each case, combined plots for different shape factors at certain 
speed/depth combination are shown. 
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Figure 6. Typical performance assessment of a submarine. 
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B. PERISCOPE SUBMERGENCE CRITERION RESULTS 



The operability indices for the periscope submergence criterion for the limited 
diameter and limited length cases in different speeds, operating depths and shape factors 
are shown in Tables I and II, and in Figure 7. Typical polar plots are shown in Figures 8 
through 30. Based on these results, the following conclusions can be drawn: 

1. Regardless of the variation of the parameters in the criterion, head seas appear 
to result in a larger number of expected criterion violations than following seas. This is 
valid for both limited length and limited diameter cases. 

2. Changes in the shape factors appear to have greater effects at smaller speeds 
and larger operating depths. At smaller depths, smaller shape factors yield slightly higher 
operability indices, while at larger depths and larger speeds, the reverse may occur. 

3. An optimum shape factor for a certain operating depth, which minimizes the 
expected number of periscope submergence events can be found, and this appears to be a 
weak function of speed. 

4. Similar values of the operability index may result in very different response 
characteristics. This is because the shape shape of the polar plots may be quite different, 
even though their areas are the same. 
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LIMITED DIAMETER 






U=3 


U=5 


U=8 


U=11 




n=2 


0.7228 


0.6614 


0.6270 


0.6154 


h=1.5D 


n=3 


0.6244 


0.6102 


0.5974 


0.5903 




n=4 


0.6192 


0.6080 


0.5935 


0.5840 




n=2 


0.7145 


0.7351 


0.6239 


0.6053 


h=2D 


n=3 


0.6772 


0.6332 


0.6163 


0.5898 




n=4 


0.6003 


0.5914 


0.5823 


0.5806 




n=2 


0.7050 


0.6309 


0.5941 


0.6109 


h=2.5D 


n=3 


0.6382 


0.6313 


0.6947 


0.6165 




n=4 


0.6897 


0.7379 


0.6578 


0.6035 



Table I : Operability indices for periscope submergence criterion for the limited diameter 
case. 







LIMITED LENGTH 






U=3 


U=5 


U=8 


U=11 




n=2 


0.6313 


0.6153 


0.6027 


0.6037 


h=1.5D 


n=3 


0.6244 


0.6102 


0.5974 


0.5903 




n=4 


0.6264 


0.6102 


0.5955 


0.5897 




n=2 


0.7717 


0.7268 


0.6320 


0.6219 


h=2D 


o=3 


0.6772 


0.6332 


0.6163 


0.5898 




n=4 


0.6266 


0.6102 


0.5973 


0.6009 




o=2 


0.7703 


0.6640 


0.6097 


0.5852 


h=2.5D 


o=3 


0.6382 


0.6313 


0.6947 


0.6165 




o=4 


0.6863 


0.7349 


0.6513 


0.6037 



Table II : Operability indices for periscope submergence criterion for the limited length 
case. 
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Oper. Index Oper. Index Oper. Index 



LIMITED DIAMETER (H=1.5D) 




■n=2 

n=3 

n=4 



LIMITED DIAMETER (H=2D) 




■n=2 

■n=3 

n=4 



Speed (Knots) 

LIMITED DIAMETER (H=2.5D) 




-n=2 

-n=3 

-n=4 



LIMITED LENGTH (H=1.5D) 
0.8r 



x 

o 

?0.7 



q.0-6 

0 



0.5 



>n=2 

■n=3 

n=4 



4 6 8 10 

Speed (Knots) 
LIMITED LENGTH (H=2D) 




■n=2 

■n=3 

n=4 



LIMITED LENGTH (H=2.5D) 




■n=2 

■n=3 

n=4 



Figure 7. 01 vs. submarine speed plots for periscope submergence criterion. 
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— n=2 

— n=3 

— n=4 
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Figure 8. Sea state-polar plot, showing the effect of shape factors in limited diameter 
case, for U = 3 Knots, h=1.5D. 
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Figure 9. Sea state-polar plot, showing the effect of shape factors in limited length case, 
for U=3 Knots, h=1.5Z>. 
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Figure 10. Sea state-polar plot, showing the effect of shape factors in limited diameter 
case, for U=5 Knots, h=1.5Z). 
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Figure 1 1. Sea state-polar plot, showing the effect of shape factors in limited length case, 
for U=5 Knots, h=1.5D. 
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Figure 12. Sea state-polar plot, showing the effect of shape factors in limited diameter 
case, for U=8 Knots, h=1.5Z). 
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Figure 13. Sea state-polar plot, showing the effect of shape factors in limited length case, 
for U=8 Knots, h=1.5ZX 
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Figure 14. Sea state-polar plot, showing the effect of shape factors in limited diameter 
case, for U=1 1 Knots, h=1.5£). 
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Figure 15. Sea state-polar plot, showing the effect of shape factors in limited length case, 
for U=ll Knots, 1.5D. 
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Figure 16. Sea state-polar plot, showing the effect of shape factors in limited diameter 
case, for U=3 Knots, h=2 D. 
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Figure 17. Sea state-polar plot, showing the effect of shape factors in limited length case, 
for U=3 Knots, h=2 D. 
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Figure 18. Sea state-polar plot, showing the effect of shape factors in limited diameter 
case, for U=5 Knots, h=2ZX 
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Figure 19. Sea state-polar plot, showing the effect of shape factors in limited length case, 
for U=5 Knots, h=2 D. 
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Figure 20. Sea state-polar plot, showing the effect of shape factors in limited diameter 
case, for U=8 Knots, h=2ZX 
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Figure 21. Sea state-polar plot, showing the effect of shape factors in limited length case, 
for U=8 Knots, h=2£>. 
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Figure 22. Sea state-polar plot, showing the effect of shape factors in limited diameter 
case, for U=1 1 Knots, h=2D. 
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Figure 23. Sea state-polar plot, showing the effect of shape factors in limited length case, 
for U=1 1 Knots, h=2D. 
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Figure 24. Sea state-polar plot, showing the effect of shape factors in limited diameter 
case, for U=3 Knots, h=2.5£>. 




270 



Figure 25. Sea state-polar plot, showing the effect of shape factors in limited length case, 
for U=3 Knots, h=2.5D. 
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Figure 26. Sea state-polar plot, showing the effect of shape factors in limited diameter 
case, for U=5 Knots, h=2.5£>. 
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Figure 27. Sea state-polar plot, showing the effect of shape factors in limited length case, 
for U=5 Knots, h=2.5D. 
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Figure 28. Sea state-polar plot, showing the effect of shape factors in limited diameter 
case, for U=8 Knots, h=2.5Z>. 




270 



Figure 29. Sea state-polar plot, showing the effect of shape factors in limited length case, 
for U=8 Knots, h=2.5D. 
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Figure 30. Sea state-polar plot, showing the effect of shape factors in limited diameter 
case, for U=ll Knots, h=2.5D. 
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Figure 31. Sea state-polar plot, showing the effect of shape factors in limited length case, 
for U=ll Knots, h=2.5D. 
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C. SAIL BROACHING CRITERION RESULTS 



The operability indices for the sail broaching criterion for the limited diameter and 
limited length cases in different speeds, operating depths and shape factors are shown in 
Tables in and IV, and in Figure 32. Typical polar plots are shown in Figures 33 through 
56. Based on these results, the following conclusions can be drawn: 

1. At smaller speeds, higher shape factors appear to result in smaller number of 
expected criterion violations. Smaller shape factors generally yield smaller indices. 

2. The operability index tends to increase with increasing depth, and in general, it 
is a weak function of speed for all shape factors. 

3. The operability index does not appear to depend on sea direction consistently. 
At certain directions, the operability index decreases significantly for various shape factor 
and speed/depth combinations. 



35 







LIMITED C 


HAMETER 






U=3 


U=5 


U=8 


U=11 




n=2 


0.6101 


0.6355 


0.7321 


0.5423 


h=1.5D 


n=3 


0.7123 


0.7138 


0.7432 


0.7013 




n=4 


0.8405 


0.8098 


0.8029 


0.7388 




n=2 


0.9091 


0.8768 


0.8969 


0.7627 


h=2D 


n=3 


0.9135 


0.9698 


0.9635 


0.9061 




n=4 


0.9235 


0.9356 


0.9333 


0.9816 




n=2 


0.8466 


0.7435 


0.8083 


0.7007 


h=2.5D 


n=3 


1 .0000 


1.0000 


0.7813 


0.8746 




n=4 


0.9777 


0.8661 


0.9421 


0.8258 



Table HI : Operability indices for the sail broaching criterion for the limited diameter case. 







LIMITED LENGTH 






U=3 


U=5 


U=8 


U=11 




n=2 


0.6173 


0.6726 


0.6942 


0.6070 


h=1.5D 


n=3 


0.7123 


0.7138 


0.7432 


0.7013 




n=4 


0.7221 


0.7314 


0.7452 


0.7717 




n=2 


0.8744 


0.9005 


0.9006 


0.8634 


h=2D 


n=3 


0.9135 


0.9698 


0.9635 


0.9061 




n=4 


0.9068 


0.9333 


0.9608 


0.9479 




n=2 


0.6925 


0.6920 


0.7799 


0.7395 


h=2.5D 


n=3 


1 .0000 


1 .0000 


0.7813 


0.8746 




n=4 


0.9761 


0.8963 


0.9408 


0.8310 



Table IV : Operability indices for sail broaching criterion for the limited length case. 
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Oper. Index Oper. Index Oper. Index 



LIMITED DIAMETER (H=1.5D) 




LIMITED LENGTH (H=1.5D) 




4 6 8 10 



Speed (Knots) 

LIMITED DIAMETER (H=2D) 




Speed (Knots) 
LIMITED LENGTH (H=2D) 




0 0.6 



4 6 8 10 



Speed (Knots) 

LIMITED DIAMETER (H=2.5D) 




Speed (Knots) 
LIMITED LENGTH (H=2.5D) 




Speed (Knots) 



Speed (Knots) 



Figure 32. 01 vs. submarine speed plots for sail broaching criterion. 
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Figure 33. Sea state-polar plot, showing the effect of shape factors in limited diameter 
case, for U=3 Knots, h=1.5Z). 
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Figure 34. Sea state-polar plot, showing the effect of shape factors in limited length case, 
for U=3 Knots, h=1.5D. 
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Figure 35. Sea state-polar plot, showing the effect of shape factors in limited diameter 
case, for U=5 Knots, h=1.5Z). 
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Figure 36. Sea state-polar plot, showing the effect of shape factors in limited length case, 
for U=5 Knots, h=1.5Z>. 
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Figure 37. Sea state- polar plot, showing the effect of shape factors in limited diameter 
case, for U=8 Knots, h=1.5ZX 
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Figure 38. Sea state-polar plot, showing the effect of shape factors in limited length case, 
for U=8 Knots, h=1.5D. 
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Figure 39. Sea state-polar plot, showing the effect of shape factors in limited diameter 
case, for U=1 1 Knots, h=1.5D. 
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Figure 40. Sea state-polar plot, showing the effect of shape factors in limited length case, 
for U=1 1 Knots, h=1.5Z>. 
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Figure 41. Sea state-polar plot, showing the effect of shape factors in limited diameter 
case, for U=3 Knots, h=2 D. 
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Figure 42. Sea state-polar plot, showing the effect of shape factors in limited length case, 
for U=3 Knots, h=2£>. 
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Figure 43. Sea state-polar plot, showing the effect of shape factors in limited diameter 
case, for U=5 Knots, h=2 D. 
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Figure 44. Sea state-polar plot, showing the effect of shape factors in limited length case, 
for U=5 Knots, h=2D. 
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Figure 45. Sea state-polar plot, showing the effect of shape factors in limited diameter 
case, for U=8 Knots, h=2D. 




Figure 46. Sea state-polar plot, showing the effect of shape factors in limited length case, 
for U=8 Knots, h=2 D. 
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Figure 47. Sea state-polar plot, showing the effect of shape factors in limited diameter 
case, for U=1 1 Knots, h=2 D. 
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Figure 48. Sea state-polar plot, showing the effect of shape factors in limited length case, 
for U=ll Knots, h=2D. 
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Figure 49. Sea state-polar plot, showing the effect of shape factors in limited diameter 
case, for U=3 Knots, h=2.5Z>. 
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Figure 50. Sea state-polar plot, showing the effect of shape factors in limited length case, 
for U=3 Knots, h=2.5ZX 
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Figure 51. Sea state-polar plot, showing the effect of shape factors in limited diameter 
case, for U=5 Knots, h=2.5Z). 
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Figure 52. Sea state-polar plot, showing the effect of shape factors in limited length case, 
for U=5 Knots, h=2.5D. 
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Figure 53. Sea state- polar plot, showing the effect of shape factors in limited diameter 
case, for U=8 Knots, h=2.5 D, 
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Figure 54. Sea state-polar plot, showing the effect of shape factors in limited length ease, 
for U=8 Knots, h=2.5D. 
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Figure 55. Sea state-polar plot, showing the effect of shape factors in limited diameter 
case, for U=1 1 Knots, h=2.5D. 
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Figure 56. Sea state-polar plot, showing the effect of shape factors in limited length case, 
for U=ll Knots, h=2.5Z>. 
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D. COMBINED CRITERIA RESULTS 



The operability indices for the combined criteria for the limited diameter and 
limited length cases in different speeds, operating depths and shape factors are shown in 
Tables V and VI, and in Figure 57. Typical polar plots are shown in Figures 58 through 
81. Based on these results, the following conclusions can be drawn: 

1. At smaller operating depths, varying the shape factors appear to yield a slight 
change in the number of expected criterion violations. At higher depths, smaller shape 
factors result in smaller operability indices. 

2. Velocity has a little effect on the operability indices at smaller operating depths. 
As the depth increases, the indices appear to decrease with increasing speed. 
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LIMITED C 


IIAMETER 






U=3 


U=5 


U=8 


U=11 




n=2 


0.4259 


0.4259 


0.4297 


0.3085 


h=1.5D 


n=3 


0.4694 


0.4689 


0.5067 


0.4492 




n=4 


0.5474 


0.5107 


0.4875 


0.4976 




n=2 


0.6236 


0.6423 


0.5307 


0.3861 


h=2D 


n=3 


0.5907 


0.6030 


0.5798 


0.4959 




n=4 


0.5238 


0.5270 


0.5179 


0.5622 




n=2 


0.5516 


0.3744 


0.4024 


0.3116 


h=2.5D 


n=3 


0.6382 


0.6040 


0.4760 


0.491 1 




n=4 


0.6674 


0.6342 


0.5999 


0.4293 



Table V : Operability indices for combined criteria for the limited diameter case. 





LIMITED LENGTH 


U=3 


11=5 


11=8 


11=11 




n=2 


0.4147 


0.4295 


0.4396 


0.3974 


h=1.5D 


n=3 


0.4694 


0.4689 


0.5067 


0.4492 




n=4 


0.4762 


0.4772 


0.4651 


0.4638 




n=2 


0.6461 


0.6408 


0.5286 


0.4994 


h=2D 


n=3 


0.5907 


0.6030 


0.5798 


0.4959 




n=4 


0.5334 


0.5435 


0.5676 


0.5488 




n=2 


0.4628 


0.3560 


0.3896 


0.3247 


h=2.5D 


n=3 


0.6382 


0.6313 


0.4760 


0.4911 




n=4 


0.6624 


0.6312 


0.5921 


0.4347 



Table VI : Operability indices for combined criteria for the limited length case. 
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Oper. Index Oper. Index Oper. Index 



LIMITED DIAMETER (H=1.5D) 

1l ; 



LIMITED LENGTH (H=1.5D) 



0.8 

0.6 

0.4 



4 6 8 10 

Speed (Knots) 

LIMITED DIAMETER (H=2D) 




^n=2 

n=3 

n=4 




n=2 

n=3 

n=4 



Speed (Knots) 

LIMITED DIAMETER (H=2.5D) 




-n=2 

-n=3 

a n=4 




*n=2 

n=3 

n=4 



6 8 10 
Speed (Knots) 

LIMITED LENGTH (H=2D) 




^n=2 

n=3 

n=4 



4 6 8 10 

Speed (Knots) 

LIMITED LENGTH (H=2.5D) 




*n=2 

■n=3 

n=4 



Speed (Knots) 



Speed (Knots) 



Figure 57. 01 vs. submarine speed plots for combined criterion. 
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Figure 58. Sea state-polar plot, showing the effect of shape factors in limited diameter 
case, for U=3 Knots, h=1.5D. 
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Figure 59. Sea state-polar plot, showing the effect of shape factors in limited length case, 
for U=3 Knots, h=1.5Z). 
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Figure 60. Sea state-polar plot, showing the effect of shape factors in limited diameter 
case, for U=5 Knots, h=1.5Z>. 
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Figure 61. Sea state-polar plot, showing the effect of shape factors in limited length case, 
for U=5 Knots, h=1.5£>. 
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Figure 62. Sea state-polar plot, showing the effect of shape factors in limited diameter 
case, for U=8 Knots, h=1.5Z). 




Figure 63. Sea state-polar plot, showing the effect of shape factors in limited length case, 
for U=8 Knots, h=1.5£>. 
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Figure 64. Sea state-polar plot, showing the effect of shape factors in limited diameter 
case, for U=1 1 Knots, h=1.5D. 
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Figure 65. Sea state-polar plot, showing the effect of shape factors in limited length case, 
for U= 11 Knots, h=1.5ZX 
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Figure 66. Sea state-polar plot, showing the effect of shape factors in limited diameter 
case, for U=3 Knots, h=2 D. 
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Figure 67. Sea state-polar plot, showing the effect of shape factors in limited length case, 
for U=3 Knots, h=2 D. 
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Figure 68. Sea state-polar plot, showing the effect of shape factors in limited diameter 
case, for U=5 Knots, h=2 D. 
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Figure 69. Sea state-polar plot, showing the effect of shape factors in limited length case, 
for U=5 Knots, h=2D. 
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Figure 70. Sea state-polar plot, showing the effect of shape factors in limited diameter 
case, for U=8 Knots, h=2 D. 
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Figure 71. Sea state-polar plot, showing the effect of shape factors in limited length case, 
for U=8 Knots, h=2 D. 
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Figure 72. Sea state-polar plot, showing the effect of shape factors in limited diameter 
case, for U=1 1 Knots, h=2£). 
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Figure 73. Sea state-polar plot, showing the effect of shape factors in limited length case, 
for U=ll Knots, h=2£>. 
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Figure 74. Sea state-polar plot, showing the effect of shape factors in limited diameter 
case, for U=3 Knots, h=2.5 D. 




270 



Figure 75. Sea state-polar plot, showing the effect of shape factors in limited length case, 
for U=3 Knots, h=2.5D. 
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Figure 76. Sea state-polar plot, showing the effect of shape factors in limited diameter 
case, for U=5 Knots, h=2.5£>. 
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Figure 77. Sea state-polar plot, showing the effect of shape factors in limited length case, 
for U=5 Knots, h=2.5D. 
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Figure 78. Sea state-polar plot, showing the effect of shape factors in limited diameter 
case, for U=8 Knots, h=2.5£>. 
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Figure 79. Sea state-polar plot, showing the effect of shape factors in limited length case, 
for U=8 Knots, h=2.5£>. 
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Figure 80. Sea state-polar plot, showing the effect of shape factors in limited diameter 
case, for U=1 1 Knots, h=2.5Z>. 
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Figure 81. Sea state-polar plot, showing the effect of shape factors in limited length case, 
for U=1 1 Knots, h=2.5D. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

In this research, the effect of geometric hull parameters, such as the length, the 
diameter and the shape factors of the hull, on the vertical plane response of submersible 
vehicles in the proximity of a free surface in deep water have been evaluated. Two cases, 
namely limited diameter and limited length, are considered, where in the first case overall 
length and in the second case the diameter is kept constant. The responses are evaluated 
by using a potential flow, strip theory solver for each case and for certain shape factors 
and different speed/operating depth combinations. Periscope submergence, sail broaching 
and the combination of both criteria are considered, in order to plot the operability 
envelope of the vehicle for different sea states and sea directions, using the above 
parameters. An operability index is calculated to quantify the operability of the vehicle. 
The main conclusions which are drawn from the study are shown below: 

1 . For the periscope submergence criterion, regardless of the variation of the parameters, 
head seas appear to result in a larger number of expected criterion violations than the 
following seas. Changes in the shape factors appear to have greater effects at smaller 
speeds and larger depths. An optimum shape factor for a certain operating depth, 
which minimizes the expected number of periscope submergence events can be found, 
and this appears to be a weak function of speed. Since the shapes of the operability 
envelopes can be quite different, similar values of the operability index may result in 
very different response characteristics. 

2. For the sail broaching criterion, smaller shape factors generally yield smaller indexes. 
The operability index tends to increase as the depth increases, and in general, it is a 
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weak function of speed for all shape factors. The operability index does not appear to 
depend on sea direction consistently. At certain directions, the index decreases 
significantly for various shape factors and speed/depth combinations. 

3. For the combined criterion, at smaller operating depths, changing the shape factors 
appears to yield a slight change in the number of expected criterion violations. At 
higher depths, smaller shape factors result in smaller operability indexes. Velocity has a 
little effect on the operability indexes at smaller operating depths. As the depth 
increases, the indexes appear to decrease with increasing speed. 

B. RECOMENDATIONS 

For further research on near surface response of submersible vehicles, the 

following studies are recommended: 

1 . In higher sea states for periscope submergence, even though the criterion is not 
violated, the average wave height may exceed the exposed periscope length, since the 
motion point appears to move more in phase with the incoming waves at higher sea 
states. This may cause difficulties in the operations, because the periscope moves in 
phase with the waves and the operator’s visual horizon may be very small. Such 
situations should be analyzed with proper simulation studies. 

2. Evaluating the effects of second order wave forces and motions on vehicle response. 
Even though these motions are slowly varying and can be controlled to a certain 
extent, they may alter both the values of the operability indices and the shape of the 
corresponding polar plots. 
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APPENDIX 



PROGRAMS 

In this Appendix, two Matlab programs for the limited diameter and the limited length 
cases and one Fortran program are presented. The purpose of the Matlab programs, as was 
mentioned in Chapter II, is to calculate the overall length for the limited diameter case, and 
the diameter for the limited length case. The programs take the shape factor values, n a and n f , 
as input and calculate the prismatic coefficients in both cases, the length in the limited 
diameter case, and the diameter in the limited length case. The Fortran program, DINIT.FOR, 
is used to create the necessary input file, SHIPMO.IN, for the strip theory seakeeping 
prediction program, SHIPMO.FOR [Ref. 2].The program reads the data written by the 
Matlab programs, and generates the input file, SHIPMO.IN. 
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% LIMDIA.M Limited diameter 

% 

% This program takes the shape factors, na and nf, 

% as input, and calculates the corresponding 
% length for the limited diameter case. 

global na nf 

d=30; 

la=3 . 6*d; 

lf=2 . 4*d; 

v=217337 . 73/ 

na=input (' Enter Na 

nf=input (' Enter Nf :'); 

xf=0 : . 3 : 1 f ; 

xa=0 : ,3:1a; 

cpa =qua d ( ' f uncpa ',0,1, .001) ; 
cpf=quad ( ' funcpf ',0,1, . 001) ; 
l=d* (4*v/pi/d / '3-3 . 6* cpa -2 . 4*cpf+6) ; 
bpl=l ; 

f id=fopen ( ' LIM2DIN' ,'w'); 

fprintf (fid, , %10.2f%10.2f%10.2f%10.2f\n', bpl, d, na , nf ) ; 
fclose (fid) ; 
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% LIMLEN.M Limited length 

% 

% This program takes the shape factors, na and nf, 

% as input, and calculates the corresponding 
% diameter for the limited length case. 

global na nf 
1=360; 
la=108; 
lf-72; 

v=217337 .73; 

na=input (' Enter Na :'); 

nf=input (' Enter Nf :'); 

xf=0 : . 3 : If ; 
xa=0 : . 3 : la ; 

cpa=quad ( 1 funcpa ',0,1); 
cpf=quad ( ' funcpf ',0,1); 
c (1) =3 . 6*cpa+2 . 4*cpf-6; 
c (2) =1; 
c (3) =0; 

c(4) = (-4*v/pi) ; 
dr=roots (c) ; 
d=dr ( 2 ) ; 

la=3 . 6*d; 
lf=2 . 4*d; 

xf=0 : . 3 : 1 f ; 
xa=0 : . 3 : la ; 

cpa=quad ( ' funcpa ',0,1); 
cpf=quad ( ' funcpf ',0,1); 
c (1) =3 . 6*cpa+2 . 4*cpf-6; 
c (2) =1; 
c (3) =0; 

c (4) = (-4*v/pi) ; 
dr=roots (c) ; 
d=dr (2) ; 
bpl=l; 

f id=f open ( ' LIM2DIN ' , ' w ' ) ; 

fprintf (fid, '%10.2f%10.2f%10.2f%10.2f\n' , bpl, d, na, nf ) 
fclose (fid) ; 
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oooooooo 



DINIT. FOR 



PROGRAM INITIALIZATION 

This program generates the SHIPMO.INI file as an input file 
for SHIPMO. FOR program. The variables are defined in the 
program manual. 



REAL BPL, LA, LF, PMB, D, NA, NF, XBKF, XBKA, BKWTH 
REAL RO, GRAV, GNU, CDISPL, DEPTH, XAXIS (21) 

REAL BILGRD, YAXIS (21, 15) , ZAXIS (21 , 15) 

REAL ZCG, RADGRO, CXCG, RADGYR, WANGI , WANGA, DWANG 
REAL XZNERT, XAC, YAC, ZAC, WA, SWL, BWL, DELWL 
REAL VMIN, VMAX, DELV, PARAM, FINISH, PI, NT, IS 
REAL XF (21) , XA(21) , YF (21) , YA(21) , CG 
INTEGER IA, IB, IC, ID, IE, IF, IG, IH, II, IJ, IK, IL, IN 
INTEGER NS, NPAC, I, ICALRD, ILID,NTO,NTI 
CHARACTER LABEL* 80 



IA=0 
IB=0 
IC=1 
ID=2 
IE=0 
IF=0 
IG=0 
IH=1 
11=1 
IJ=0 
IK=0 
IL=0 
IN=0 
NS=2 1 
NPAC= 1 
BKWTH=1 . 0 
XBKA=- 10.0 
XBKF=10 . 0 
RO=l .9905 
GRAV=32 .1740 
GNU=1 . 26E-05 
CDISPL=0 . 0 
DEPTH=0 . 0 
BILGRD=0. 0 
ILID=0 
RADGRO=0 . 0 
CXCG=9999 . 0 
RADGYR=0 . 0 
XZNERT=0 . 0 
ICALRD=5 
FINISH=0 . 0 
PI=3. 1415927 
NTO=15 
NTI=1 
WA=1 . 0 
SWL=0 . 2 
BWL=2 . 0 
DELWL=0 . 1 
WANGI=0 . 0 
WANGA=345.0 
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DWANG=15. 0 

BPL=360 . 0 
D=30 . 0 
NA=3 . 0 
NF=3 . 0 

VMIN=5* 1 . 689 
VMAX=5* 1 . 689 
DELV=0 . 0 



Z=2 . 5*D 



OPEN (10, FILE=' SHIPMO. IN' ) 
WRITE (*,*) 'ENTER LABEL' 
READ (*, ' (A) ' ) LABEL 

OPEN (100, FILE= ' LIM2DIN ' ) 
READ (100, 200) BPL, D, NA, NF 



LA=3 . 6*D 
LF=2 . 4*D 
PMB=BPL-6*D 



ZCG= (-1) *Z- (D/2) -3.0 
CG= (D/2) -3 . 0 
YAOO.O 

XAC= (BPL/ 2 . ) - (BPL/3 . ) 

ZAC= (-1 ) *Z+D 

WRITE (*,*) 'ENTER PARAM' 

READ *, PARAM 1 

PARAM=PARAM 1 *3.28 
WRITE (10,*) LABEL 

WRITE (10, 101) IA, IB, IC, ID, IE, IF, IG, IH, II, IJ, IK, IL, IN,NS,NPAC 
WRITE (10, 1000) BPL, D,NA,NF, PARAM, Z,CG 
WRITE (10, 102) BPL, RO, GRAV, GNU, CDISPL, DEPTH 
WRITE (10, 103) XBKF, XBKA, BKWTH 

DO 11 IS=1, 10 

XF (IS) =LF* (10-IS) / 9 
XAXIS ( IS) = ( PMB/ 2 ) +XF ( IS ) 

YF ( IS) = (D/2 ) * (1- (XF(IS) /LF) **NF) ** (1/NF) 

DO 12 NT=1 , 7 

Z AX IS (IS, NT ) = ( - 1 ) *Z- (D/2) -YF(IS) *SIN( (PI/2) * (8-NT) /7.) 
YAXIS (IS, NT) =YF (IS) *COS ( (PI/2) * (8-NT) /7 . ) 

12 CONTINUE 

ZAXIS (IS, 8) = (-1) *Z- (D/2) 

YAXIS (IS, 8) =YF(IS) 

DO 13 NT=9, 15 

ZAXIS (IS, NT) = (-1) *Z- (D/2) +YF(IS) *SIN( (PI/2) * (NT-8) /7. ) 
YAXIS ( IS , NT ) =YF ( IS ) *COS ( (PI/2) * (NT-8) /7 . ) 

13 CONTINUE 
11 CONTINUE 

ZAXIS (10, 15) =0 . 0 

YAXIS (10, 15) =YF (10) *COS ( (PI/2) * (6./7. ) ) 

XAXIS (11) =0.0 
DO 21 NT=1 , 7 

ZAXIS (11, NT) = (-1) *Z- (D/2) - (D/2) *SIN( (PI/2) * (8-NT) /7 . ) 
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YAXIS ( 11, NT) = (D/2 ) *COS ( (PI/2) * (8-NT) /7.) 

21 CONTINUE 

ZAXIS (11, 8)= (-1) *Z- (D/2) 

YAXIS (11, 8)=D/2 
DO 31 NT=9, 14 

ZAXIS (11,NT)=(-1)*Z- (D/2) + (D/2)*SIN( (PI/2)* (NT- 
YAXIS (11, NT) = (D/2) *COS ( (PI/2)* (NT-8)/7.) 

31 CONTINUE 

ZAXIS (11, 15)=0.0 

YAXIS (11, 15) = (D/2) *COS ( (PI/2) * (6./7. ) ) 

DO 41 IS=12, 21 

XA (IS) =LA* (IS-12 . ) / 9 . 

XAXIS ( IS ) = (-1 ) * (PMB/2) -XA(IS) 

YA ( IS) = (D/2 ) * (1- (XA(IS) /LA) **NA) 

DO 42 NT=1, 7 

ZAXIS ( IS , NT) = (—1) *Z— (D/2) -YA(IS) *SIN( (PI/2 
YAXIS (IS, NT) =YA( IS) *COS ( (PI/2)* (8-NT)/7.) 
42 CONTINUE 

ZAXIS (IS, 8) = (-1) *Z- (D/2) 

YAXIS (IS, 8)=YA(IS) 

DO 43 NT=9, 15 

ZAXIS (IS, NT) = (-1) *Z- (D/2) +YA(IS) *SIN( (PI/2 
YAXIS (IS, NT)=YA(IS) *COS ( (PI/2)* (NT-8)/7.) 
CONTINUE 



NTI, XAXIS (1) , BILGRD, ILID 
YAXIS (1, 1) , ZAXIS (1, 1) 



43 

41 CONTINUE 

WRITE (10, 104) 

WRITE (10, 105) 

DO 51 IS=2 , 2 0 

WRITE (10,104) NTO,XAXIS (IS) , BILGRD, ILID 
DO 52 NT=1, 15 

WRITE (10, 105) YAXIS (IS, NT) , ZAXIS (IS, NT) 
52 CONTINUE 



51 CONTINUE 

WRITE (10, 104) 
WRITE (10, 105) 
WRITE (10, 106) 
WRITE (10, 107) 
WRITE (10, 108) 
WRITE (10, 109) 
WRITE (10, 110) 
WRITE (10, 111) 
WRITE (10, 112) 



NTI, XAXIS (21) , BILGRD, ILID 
YAXIS (21,1), ZAXIS (21, 1) 

ZCG, RADGRO 

CXCG, RADGYR, XZNERT 

ICALRD, XAC, YAC, ZAC 

WA, SWL, BWL, DELWL, VMIN, VMAX, DELV 

WANG I , WANGA, DWANG 

PARAM 

FINISH 



101FORMAT (1515) 

102 FORMAT (3F10 . 4, 2E15 . 5, F10.4) 
103 FORMAT ( 3 FI 0 . 4 ) 

104 FORMAT (15, 2 FI 0. 4, 15) 

1 05 FORMAT (2F10.4) 

10 6 FORMAT (2F10. 4) 

107 FORMAT (2F10 . 4, F15 . 3) 
108FORMAT (12, F8.4,2F10.4) 
109FORMAT (8F10. 4) 

110 FORMAT ( 3 FI 0 . 4 ) 

111 FORMAT ( FI 0 . 4 ) 

112 FORMAT ( F4 . 1 ) 

200 FORMAT ( 4 FI 0 . 2 ) 

1000 FORMAT ( 7 FI 0 . 2 ) 

END 



8)/7.) 



) * (8-NT) /I .) 



) * (NT-8) /I . ) 
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